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We propose to apply spin noise spectroscopy (SNS) to detect many-body localization (MBL) in 
disordered spin systems. The SNS methods are relatively non-invasive technique to probe sponta¬ 
neous spin fluctuations. We here show that the spin noise signals obtained by cross-correlation SNS 
with two probe beams can be used to separate the MBL phase from a noninteracting Anderson 
localized phase and a delocalized (diffusive) phase in the studied models for which we numerically 
calculate real time spin noise signals and their power spectra. For the archetypical case of the 
disordered XXZ spin chain we also develop a simple phenomenological model. 


The fate of Anderson localization in the presence 
of inter-particle interactions in a disordered quantum 
medium is an exciting frontier in condensed matter 
physics [THIS]- It is well known that all states of a one¬ 
dimensional (ID) disordered chain of noninteracting par¬ 
ticles are Anderson localized (AL) for any amount of dis¬ 
order [20l[2T]. Thus, the AL state is a perfect insulator as 
long as the particles are not coupled to other degrees of 
freedom. In the presence of inter-particle interactions, a 
dynamical transition from a delocalized (diffusive) phase 
to a many-body localized (MBL) phase has been pre¬ 
dicted in disordered quantum media when the strength 
of (quenched) randomness is increased [4]-[9]. 

The MBL state is also a perfect insulator, but it has 
different dynamical properties compared to the AL state 
of noninteracting particles [TOHIIIIS]- For example, en¬ 
tanglement entropy in the MBL phase shows a slow log¬ 
arithmic growth following a global quench in an isolated 
system [TOHE]- However, the entanglement entropy is 
very difficult to measure experimentally, and electrons or 
spins in conventional solid-state systems are coupled to 
an environment such as a phonon bath. Recent studies 
show that signatures of the MBL phase can survive in the 
presence of weak coupling to a thermalizing environment 
[22] [23]. In particular, the spectral functions of local op¬ 
erators can be used to identify the MBL phase in the 
presence of weak dissipation [22l [23] . 

A new experimental approach based on a modified non¬ 
local spin-echo protocol along with a double electron- 
electron resonance technique in electron spin resonance 
has been proposed to distinguish the MBL phase from 
a noninteracting AL phase and a delocalized phase at 
infinite temperature [24]. The proposed approach can 
probe interaction effects, thus is able to separate the 
MBL phase from the AL phase. The method does in¬ 
volve optical pumping or polarization of local spins by 
external pulse fields, which can in principle lead to un¬ 
wanted local heating and excitations. Here we propose a 
relatively non-invasive method based on spin noise spec¬ 
troscopy (SNS) to distinguish the MBL phase from the 
AL phase and the delocalized phase in disordered spin 
systems. 

The optical SNS method has been developed recently 
as an alternative to conventional perturbation-based 
(pump-probe) techniques for measuring dynamical spin 



FIG. 1 . (a) Schematic of the two-beam spin noise spec¬ 

troscopy set-up with two probe lasers detecting spin fluctua¬ 
tions of different spins, (b) Disorder (77) vs. energy density 
(e) phase-diagram of the disordered Heisenberg chain EH. 


properties ESI US). Intrinsic spin fluctuations of electrons 
and holes are passively detected in SNS by measuring op¬ 
tical Faraday rotation fluctuations of a linearly polarized 
probe laser beam passing through the sample [27H38] . 
SNS with a single probe laser has been useful in char¬ 
acterizing various properties (e.g., g-factors, relaxation 
rates and decoherence times) of different spin ensembles, 
such as specific alkali atoms m, itinerant electron spins 
in semiconductors [33] and localized hole spins in quan¬ 
tum dot ensembles [34]. Last year, an extension of the 
traditional SNS has been proposed and demonstrated by 
using two linearly polarized probe lasers for detecting 
inter-species spin interactions in a heterogeneous two- 
component spin ensemble interacting via binary exchange 
coupling in thermal equilibrium PIEO]. In this cross¬ 
correlation SNS method, intrinsic spin fluctuations from 
two different species are independently detected, and in¬ 
teraction effects are determined by the cross-correlation 
of these two spin noise signals. Interaction effects be¬ 
tween spins of a single species in a sample can also 
be determined by two-beam SNS when the two probe 
laser beams are spatially separated as shown in Fig. Ba) 
Elm]- Thus, we propose that the two-beam SNS mea¬ 
surements can be used to separate the MBL phase from 
an AL phase. In fact, as we show in our example, the 
responses from single-beam and two-beam SNS can be 
efficiently employed to distinguish the MBL phase from 
both the AL and delocalized phases. 

In order to demonstrate our idea we study the spin- 
noise signals of a ID random-field XXZ spin chain. 
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perhaps the best-studied “canonical” model system for 
MBL. In addition, results for the disordered transverse- 
field Ising model are given in the Suppl. Mat. [43]. The 
Hamiltonian of disordered XXZ spin chain is given by 

N 

H = + hiS!,{i) 

where the random fields hi are chosen from a uniform dis¬ 
tribution within the window [— 77 , 77 ], and 
for periodic boundary conditions. For isotropic exchange 
interactions {Jz = J±), this is a disordered Heisenberg 
chain. The spin Hamiltonian in Eq. can be mapped to 
a model of interacting spinless Jordon-Wigner fermions 
on a lattice where J±_ represents hopping between neigh¬ 
boring sites, Jz is inter-particle interaction strength and 
hi denotes a random local chemical potential [42|. For 
Jz = 0^ the Hamiltonian in Eq. represents noninter¬ 
acting fermions subject only to a random local poten¬ 
tial, and hence can be used to study the AL phase of 
noninteracting particles. This model with a non-zero Jz 
has been extensively studied recently in the context of 
MBL. A disorder versus energy density phase-diagram 
separating the delocalized and MBL phases of an iso¬ 
lated chain is shown in Fig.j^b) following Ref. |44| where 

6 = ((-H) -Emax)/(-^min -E'max) with E/max (E/min) being 

maximum (minimum) energy. 

For the Hamiltonian in Eq. we now describe the 
signals measured in the SNS set-ups. The measured re¬ 
sponse of the spin component tranverse to the random 
magnetic field in the single-beam SNS is 

cnt) = {{{sf{t),sno)})), ( 2 ) 

where is the anti-commutator and ((..)) denotes 

both (canonical) thermal and disorder averaging. Here, 
the x-component of total spin polarization at the mea¬ 
surement spot at time t is Sj{t) = where the 

sum runs over all local spin sites within the sample spot 
illuminated by the probe laser beam /. The correlation 
function Cj{t) describes the relaxation of spontaneous 
spin fluctuations at the spot of the sample probed by 
the single-beam SNS. The single-beam (local) spin-noise 
power spectrum is obtained by its Fourier transform, 

/ CX) 

dt (3) 

-00 

Let \En) denote a many-body eigenstate with eigenvalue 
H\En) = En\En). The single-beam spin-noise power 
spectrum reads, 

Pfi^) = (4 

n,p 

\{Er.\Sf\Ep)\^), (4) 

with partition function Z = ^ inverse tem¬ 

perature P = l/ksT. (..) in Eq. 1^ describes averaging 


over different disorder realizations. The cross-correlation 
signal of the two-beam SNS is 

C^rit) = {{{Sf{t),SUO)}E{Sfj{t),SnO)})). ( 5 ) 

where ^ 

through all sites which the probe laser II illuminates. 
The cross-correlation spin noise power 

/ oo 

dt ( 6 ) 

-00 

= ( I X ) 

n,p 

Re{{En\Sf\Ep){Ej,\Sfj\En))). (7) 

We calculate the eigenvalues and eigenstates of the 
Hamiltonian using exact diagonalization, and evaluate 
the spin noise signals in Eqs. |2|5| and the corresponding 
power spectra in Eqs. |4|7| In the following we present 
results using periodic boundary conditions at high tem¬ 
perature, ksT = 50J_l, averaging over 3000 disorder real¬ 
izations. However, the main results of this paper are also 
valid at moderate temperatures m- We quote Jz^hi^r] 
in the units of Jx, and fix Jx = 1 throughout the paper. 

The single-beam spin noise response in real time for 
the interacting delocalized phase at low disorder (77 = 1 ) 
and for the MBL phase at high disorder (77 = 5, at which 
every many-body eigenstate of the isolated chain is local¬ 
ized) of the random XXZ model is shown in Fig. [^a,b). 
The single-beam responses are clearly different in the two 
phases. The transverse spin component relaxes exponen¬ 
tially in the delocalized phase, while in the MBL phase 
it appears to do so algebraically, and it also oscillates. 
We show the single-beam spin noise response of the AL 
phase {Jz = 0) in Fig. |^c), which again indicates an 
algebraic relaxation of the transverse spin component. 
The relaxation in this phase can be approximated as 
sin(t/rA)/(t/rA) with some characteristic time scale 
We use high disorder (77 = 5) in Fig. |^c) to ensure that 
all states are localized even for a finite-length noninter¬ 
acting chain. We provide log scale plots in the Suppl. 
Mat. [43] for the real-time spin noise responses in Fig. 
to highlight the nature of their respective decays. 

Within the single-beam SNS measurements, the slow 
spin relaxation due to interactions in the highly disor¬ 
dered XXZ model is masked by other mechanisms of re¬ 
laxation (e.g., due to random fields and hopping) also 
present in the noninteracting case. Thus, it becomes dif¬ 
ficult to separate the MBL phase from the AL phase by 
single-beam SNS. However, interaction effects are nicely 
diagnosed by two-beam SNS when we cross correlate two 
different noise signals to exclude their self-correlations. 

The two-beam cross-correlation spin noise responses 
in real time for the three different phases are shown 
in Fig. [^d,e ,f) when two beams are just next to each 
other without overlapping (i.e., the separation between 
the centers of the two beams is two lattice spacings). 
The responses here are quite different in all three phases. 
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FIG. 2 . Single-beam and two-beam spin noise responses in 
real time in the three different phases, (a-c) Cf{t) and (d-f) 
Ccr(t) are obtained by exact diagonalization in the delocalized 
{Jz = 1,?7 = 1), the MBL (J^ — = 5) and the Anderson 

localized {Jz = 0, ?7 = 5) phase. Here = 12 and the beam 
I and II illuminate respectively spins 5,6 and 7,8. 


A large cross-correlation between the transverse compo¬ 
nents of different spins is developed in the diffusive phase, 
and it relaxes exponentially fast. The cross-correlation 
in the MBL phase is relatively smaller and it relaxes 
slowly with many oscillations. In the AL phase, the 
transverse component hardly shows any cross-correlation 
between different spins at high temperature. Notice 
that Cj(0) = 1 while C^y{ 0) ^ 0 in Fig . because 
(5'f (0)5'^(0)) = 6im/^ at high temperature j^ . We also 
perform these calculation for the relaxation of spins along 
the ^-direction, and find that it is similar in the MBL and 
AL phases even within the two-beam SNS measurements 
[43l|46]. Thus, 2 ;-direction response is not a good candi¬ 
date to separate these two phases. 

Turning to the frequency domain, the disorder aver¬ 
aged spin noise power spectra of single-beam and two- 
beam SNS are shown in Fig.[^ We find a Lorentzian line- 
shape for the single-beam noise power spectrum Pf{uj) 
[Fig-i a)] in the diffusive regime refiecting the exponen¬ 
tial spin relaxation. In the MBL phase, the shape of 
Pf{(jj) exhibits a plateau as shown in Fig. ^b). The 
overall rectangular shape of Pf{uj) in the ^L regime 
[Fig-i c)] is somewhat similar to that of the MBL phase 
but it shows much higher fiuctnations. 

Crucially, the cross-correlation spin noise power spec¬ 
tra Pcri^) of fbe two-beam SNS in Fig. d-f) differ sig¬ 
nificantly in shape and magnitude in the delocalized (dif¬ 
fusive), MBL and AL phases. The power spectrum in the 
delocalized phase [Fig. |^d)] can be perceived as the dif¬ 
ference between two equal area Lorentzians of different 
widths. Thus, the total area under the cross-correlation 
curve is zero, a signature of high-temperature behaviour. 
The line-shape of P^^{oo) in the MBL phase is very dif¬ 
ferent from that in the delocalized and AL phases. It 


FIG. 3. Single-beam and two-beam spin noise power spectra 
in the three different phases, (a-c) Pf{uj)^ (d-f) Pcr{^) are 
obtained by exact diagonalization and (g-i) Pcr{^) are calcu¬ 
lated from Eq. 10 in the delocalized {Jz — 1, = 1), the MBL 

{Jz = 1,?7 = Sjand the Anderson localized {Jz = 0,?7 = 5) 
phase. We use 7 = 0.068 ,71 = 0.18 in (g), 7 = 0.03 ,71 = 0.08 
in (h) and 7 = 0 in (i). The parameters in (a-f) are the same 
as in Fig. 


is also almost one order of magnitude smaller than in 
the delocalized phase but is clearly non-zero [Fig. [^e)]. 
However, P^A^) fbe AL phase is very small (one or¬ 
der of magnitude smaller than that in the MBL phase at 
high temperature) and featureless indicating essentially 
no correlation between spins. Note that P^i^) shows 
both negative and positive values which signify (anti- 
) correlations between spins at different frequency/time 
scales. The anti-correlations in Fig. [^d,e) at higher fre¬ 
quencies (shorter time) are induced by fast co-flips be¬ 
tween different spins in the presence of spin-exchange 
coupling. This feature is present in the two-beam cross¬ 
correlation spectra in the delocalized and MBL phases 
but is absent in the AL phase. 

Pcri^) depends on the separation between the two 
probe beams. With an increasing separation between 
the beams, P^i^) exhibits more oscillations in the delo¬ 
calized and MBL phases, and its magnitude falls rapidly 
in the MBL and AL phases as might be expected in the 
presence of spatial localization. Also the shape of P^A^) 
in the three phases remains unchanged when we explicitly 
couple the disordered XXZ chain to a weakly thermaliz- 
ing environment. Both of these features are discussed in 
the Suppl. Mat. [43]. 

In order to understand the cross-correlation noise 
power spectra obtained from numerics, we develop a sim¬ 
ple phenomenological model inspired by the experimental 
set-up of the two-beam SNS. Let us denote the respec¬ 
tive total spin polarizations of the spins illuminated by 
the probe laser beams I and // by S/ and S//. We as¬ 
sume that the spins in beam I {II) can relax (a) due to 
spin-exchange interactions with the spins in beam II (/), 
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and (b) due to interactions with other spins in the sam¬ 
ple. The spin-exchange interactions between the spins 
from the two illuminated spots conserve total spin but 
do transfer spins between the two spots via the spin co- 
flips with a rate 7 . This leads to cross-correlations of 
their fluctuations. For process (b) we define a net spin 
relaxation rate 7 ck for spins in beams a = /,//. This 
process does not conserve the total spin in the beams. 
Combining these with Larmor precession in the random 
magnetic field, we obtain the following stochastic evolu¬ 
tion equations [40] . 

^ = S„xh„-7„S„-7(S„-S5)+^„, (8) 


where 7 = //, // = /. Here we have included stochas¬ 
tic fluctuations For simplicity, we consider these as 
Gaussian white noise with zero mean. 


= S{t - i') Y {^oc/Sla + liSa/3 “ '5a^))(^) 

with i^j = x^y^z. The relation in Eq. [^ensures station¬ 
ary equal-time correlation {Sai{t)Spj{t)) oc in 

the high temperature limit m- We coarse-grain the ran¬ 
dom magnetic fields within the spot of each laser beams, 
and represent them by a single random field. Here we 
choose the random field along ^-direction following the 
Hamiltonian in Eq. i.e., = h^z where ha are two 

i.i.d. random numbers chosen from a uniform distribu¬ 
tion over [— 7 , 7 ]. Averaging over noise, we find the cross¬ 
correlator of spin polarizations in Eq.j^in compact form, 

w = iT.{z^)’ (i») 


where = 7172 - (cj ± hi){uj ± h//), k,± = 73 ( 2 cl; ± 
{hi + hii)). Here 7 / = 7 // = 71 , 72 = 71 -h 27 and 
73 = 71+7 are characteristic rates that influence broad¬ 
ening of the peaks. (..) in Eq. 10 denotes averaging over 


the random fields ha- Due to the presence of the ran¬ 
dom field throughout the sample, the phenomenological 
relaxation rates 7 , 7 / and 7 // vary with disorder strength 
7 , and we use them as fitting parameters. We plot the 
cross-correlator after averaging over the random 

fields for 77 = 1,5 in Eig. |^g,h). We find remarkably 
good agreement between the numerical results obtained 
at high temperature and the results from our phenomeno¬ 
logical model. Erom Eq. 10, P^vi^) = 0 for 7 = 0, in the 
AL phase (J^ = 0). We note, however, that our model 
assumes instantaneous transfer of polarization density 
through direct co-flips of spins, appropriate for directly 
adjacent spots. If this condition is not fulfilled, a delay 
for the propagation of polarization density over the dis¬ 
tance between the illuminated spots will have to be added 
to our simplified model. Secondly, our model works only 
at relatively high temperatures and the results obtained 
from the phenomenological model would deviate from the 
numerical results at low temperature [431147]. 

To summarize, we have calculated spin noise signals 
in the diffusive, AL and MBL regimes of the disordered 
XXZ and transverse-field Ising spin chains, for differ¬ 
ent temperatures as well as beam separations. We have 
shown how these signals can be used to identify the pres¬ 
ence of MBL. The current proposal addresses disordered 
spin systems beyond previously addressed regimes of in¬ 
finite temperature and strictly isolated systems. Eor 
an actual experiment, it is necessary to have a disor¬ 
dered spin system whose interaction energies and relax¬ 
ation time scales are within the bandwidth of the recent 
optical SNS measurements (nearly 100 GHz) [36]. In 
this regard, solid state systems, such as localized spin 
defects in solids (e.g., nitrogen-vacancy centers in dia¬ 
mond) [24l l48ll5Q] . are suitable candidates for SNS mea¬ 
surements to detect MBL. Apart from SNS, spatially re¬ 
solved spin noise from quasi-ID or 2D layers of inter¬ 
acting spins can also be measured using recently demon¬ 
strated sensitive nanoscale magnetic resonance imaging 
using nitrogen-vacancy centers in diamond [5ll|52]. 

We thank Vadim Oganesyan and Shivaji Sondhi for 
helpful discussions. 
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FIG. SI. Log scale plots of single-beam and two-beam spin 
noise responses of transverse spin component in real time in 
three different phases of disordered XXZ spin chain. (a,d) in 
log-linear scale in the delocalized phase {Jz = 1,t] = 1 ) and 
(b,c,e,f) in log-log scale in the MBL (J^ = 1,77 = 5) and the 
Anderson localized {Jz = 0,?7 = 5) phase. Here iV = 12 and 
the beam / and II illuminate respectively spins 5,6 and 7,8. 


FIG. S 2 . Single-beam and two-beam spin noise responses of 
longitudinal spin component in real time in the delocalized, 
MBL and Anderson localized phases of disordered XXZ spin 
chain. The parameters are the same as in Fig. |S1| 


We show plots in log-linear and log-log scales of real¬ 
time spin noise responses of the transverse spin com¬ 
ponent of disordered XXZ spin chain to elucidate the 
nature of their relaxation. We plot absolute values of 
Cj{t) and Ccr(t) in log scale and time t in linear scale 
in Fig. [ST]^a,d). Both the real-time responses in the de- 
localized^diffusive) phase appear to exhibit exponential 
decay in time. In Fig. SI ^b,c,e,f) we plot both absolute 
values of Cj{t) and and time t in log scale. While 

the single-beam spin noise response of the transverse spin 
component in real time shows algebraic (power-law) de¬ 
cay both in the MBL and AL phases, cross-correlation 
spin noise response shows nearly algebraic decay only in 
the MBL phase. The transverse spin component shows 
significantly less cross-correlation in the AL phase, as is 
evident in Fig. IHf)- 


In Fig. we show single-beam and two-beam spin 
noise responses of the longitudinal spin component {z- 
direction) of disordered XXZ spin chain in real time in the 
three different phases. Both spin noise responses show 
different relaxation behaviour in the delocalized and the 
localized phases. However, the relaxation of the longi¬ 
tudinal spin component is quite similar in the MBL and 
the AL phases within both single-beam and two-beam 
cross-correlation measurements. 


Next, we discuss how cross-correlation spin noise spec¬ 
tra Pcri^) of transverse spin component of the disordered 
XXZ spin chain depend on temperature and separation 
between the probe beams in the three different phases. 
For this we change temperature and separation between 
the probe beams from those in Fig. 3 of the main text. 
For comparison in Fig. |S3[d-f) we use the same temper¬ 
ature and separation between the probe beams as in the 
main text. We reduce temperature to 10Jx in Fig. S3 A- 


c), and we find that the line-shapes of P^ri^) three 

phases remains similar to that in Fig.[S^d-f) at high tem¬ 
perature 50J^. Thus, the two-beam SNS is still a good 
probe to separate the three different phases. However, 
fluctuations in P^ri^) i^ phase at low tempera¬ 

ture are increased by one order of magnitude from that at 
high temperature. In fact, Pcri^) is similar in magnitude 
both in the MBL and AL phases but still has different 
shapes. In Fig. S3^g-i) we move the beam II by one 
lattice spacing, and we find extra oscillations in P^ri^) 
compared to Fig. [S^d,e) in the delocalized and MBL 
phases. We also notice that the magnitude of P^ri^) 
is reduced by one order of magnitude in the MBL phase 
and by two orders in the AL phase compared to Fig.jS^e- 
f). We use periodic boundary conditions and averaging 
over 3000 disorder realizations in Figs. |Sl|S2|S3l except 
in Fig. |S3[h) where 10000 disorder realizations are used. 
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FIG. S3. Dependence of cross-correlation spin noise power 
spectra Pcr{^) on temperature and separation between the 
probe beams in the delocalized, MBL and Anderson localized 
phases of the disordered XXZ spin chain. Here iV = 12 and 
temperature is 10 J± in (a-c) and 50 in (d-i). The beam I 
and II illuminate respectively spins 5,6 and 7,8 in (a-f) and 
spins 5,6 and 8,9 in (g-i). 


FIG. S4. Gross-correlation spin noise power spectra Pct{^) 
in the delocalized (J^ = 1,77 = 1 ), the MBL (J^ = 1,77 = 4) 
and the Anderson localized phases {Jz = 0,?7 = 4) of the 
disordered XXZ spin chain coupled to a thermalizing bath. 
The bath coupling Jc = 0.15 in (a-c) and Jc = 0.25 in (d-f). 
Here A" = 6 , J_l = = 775 = 1, temperature is 50 J_l and 

beams I and II illuminate respectively spins 2,3 and 4,5. 


We check the nature of cross-correlation spin noise 
power spectra of disordered XXZ spin chain by explic¬ 
itly coupling the chain (system) to a thermalizing bath. 
We model the bath here by a weakly disordered XXZ 
spin chain in the diffusive phase (nonintegrable). The 
Hamiltonians of the bath and system-bath coupling are 
respectively, 


N 

bath — ^ ^ Jb^i'^i-\-l T 1 

(SI) 

N 

i=l 

(S 2 ) 


where the random fields h\ are chosen from a uniform 
distribution within the window [— 77 ^, 775 ], and aw+i = 
a I for periodic boundary conditions. Here Jc deter¬ 
mines strength of the system-bath coupling, and the bath 
should be effective in thermalizing the system when the 
coupling matrix element is of the order of the many-body 
level spacing in the bath. The full Hamiltonian of the sys¬ 
tem plus bath is i7f = + i^bath + He where H is given 

in Eq.[^of the main text. We calculate spin noise signals 
in an eigenstate \n) of Hf at an energy En corresponding 
to that of the thermal ensemble at inverse temperature 
/3, En = (n|i7f|n) = Tr[e-^^fi7f]/Tre-^^L InFig.l^we 


show two-beam cross-correlation spin noise spectra of the 
transverse spin component in the three different phases 
for two different values of Jc. Noise spectra are shown 
after averaging over 30000 (except 10000 in Fig. S4 A,d)) 
disorder realizations. We find that the shapes of two- 
beam spin noise spectra in Fig. are similar to those 
obtained earlier without explicitly coupling the system 


to a thermalizing bath. This shows the robustness of 
two-beam spin noise signals to distinguish the three dif¬ 
ferent phases. We note from Fig.j^that the fluctuation 
in the noise spectra falls with increasing bath coupling 
as long as the coupling is smaller than the characteristic 
energy scales in the system. It signals better thermal- 
ization with increasing bath coupling without destroying 
the signatures of the three different phases in the isolated 
disordered XXZ chain. We mention that the large fluc¬ 
tuations in the noise spectra shown in Fig. are also 
due to relatively shorter chain length which we were able 
to simulate for this case. 

Finally, we present cross-correlation spin noise signals 
measured by two-beam SNS for a disordered transverse- 
field Ising model with next-nearest neighbor coupling. 
The Hamiltonian of this model is [1] 


AT-l N-2 

i=l i=l 


J2 X 


i+2 


N 

hJ2Sf(,S3) 

i=l 


where the nearest neighbor couplings Ji = J + 6Ji^ with 
all random SJi chosen from a uniform distribution within 
the window [— 77 , 77 ]; h is a constant magnetic field. For 
finite J 2 , the above model has a delocalized phase at low 
disorder and an MBL phase at high disorder [T] . The sys¬ 
tem is in the AL phase when J 2 = 0. We numerically cal¬ 
culate cross-correlation spin noise signals of the spin com¬ 
ponent along 7 /-direction. We quote J 2 ,h,r] in the units 
of J, and fix J = 1. The cross-correlation spin noise sig¬ 
nals in real time in the delocalized, MBL and AL phases 
are shown in Fig. S5[a-c) at high temperature. A large 
cross-correlation between the transverse spin components 
of separate spins is developed in the delocalized phase 



































= 0 . 3,77 = 1 


J 2 = 0 . 3,77 = 5 J 2 = 0,77 = 5 



UJ U) U) 


and it relaxes relatively fast, while the cross-correlation 
is much weaker in the MBL phase and it relaxes very 
slowly. The transverse spin component hardly shows any 
cross-correlation between different spins at high tempera¬ 
ture in the AL phase. We plot cross-correlation spin noise 
power spectra in Fig.[S^d-f) which are different in shape 
and magnitude in the three phases. Especially, the cross¬ 
correlation power spectra in the MBL and AL phases are 
very different in shape. We use periodic boundary con¬ 
ditions, averaging over 3000 (in Fig. [S^a-c)), 6000 (in 
Fig. [S^d)) and 20000 (in Fig. [S^e,f)y~disorder realiza¬ 
tions. 


FIG. S5. Two-beam spin noise responses of transverse 
spin component in the three different phases of disordered 
transverse-field Ising chain, (a-c) C^r(t) and (d-f) Pcr(^) are 
obtained numerically in the delocalized (J 2 = 0 . 8,77 = 1 ), 
the MBL (J 2 = 0 . 3,77 = 5) and the Anderson localized 
(J 2 = 0,77 = 5) phase. Here A = 12, J = 1, /i = 0.6, tem¬ 
perature is 50 J and beams I and II illuminate respectively 
spins 5,6 and 7,8. 
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